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Using a general classification of dark enegy models in four classes, we discuss the complementarity 
of cosmological observations to tackle down the physics beyond the acceleration of our universe. We 
discuss the tests distinguishing the four classes and then focus on the dynamics of the perturbations 
in the Newtonian regime. We also exhibit explicitely models that have identical predictions for a 
subset of observations. 
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I. INTRODUCTION 

The flow of new data nas allowed to nar- 

row the constraints on the standard cosmological param- 
eters. Among all the conclusions concerning our universe, 
the existence of a recent acceleration phase seems to be 
more and more settled. Even though, cosmology has a 
minimal standard model, consistent and robust with 6 or 
7 free parameters (the concordance model) , many exten- 
sions (both of the primordial physics and of the matter 
content) are still weakly constrained today. These pa- 
rameters start to be measured very accurately but it may 
turn out that our (successfull) parameterization may be 
too simple. 

The quest for the understanding of the origin of this 
acceleration is however just starting (see e.g. Refs. 0, 
H,EJE3) Different ways of attacking this problem have 
been proposed 0, El 113 using various properties to or- 
der the questions on the origin of the acceleration. In 
this article, we come back to the classification proposed 
in Refs. 0, 0] and detail the degeneracies that have 
to be taken into account before conclusions are drawn. 
This completes other attempts to define strategies to get 
insight into the physics of dark energy 1 . 

The conclusion that our universe is accelerating as- 
sumes in the first place the validity of the Copernician 
principle, that is the fact that the observable universe 
can be described on large scales by a Friedmann-Lemaitre 
spacetime so that its dynamics is characterized by a single 
function of time, the scale factor a(t). Since most data 
are related to events observed on our past light cone, 
there is an intrinsic degeneracy along this cone. In the 
Friedmann-Lemaitre models, this degeneracy is lifted by 
the symmetry assumption. In the next to simple case, 
the universe can be described by a spherically symmetric 
spacetime. For an observer seating at the center of the 
universe, the degeneracy along the past-null cone entan- 
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1 We emphasize here that the cause of the acceleration of the uni- 
verse is decoupled from the cosmological constant problem which 
is somehow assumed to be solved. One goal of the dynamical 
dark energy models is to avoid the fine tuning related to the 
cosmic coincidence problem. 



gles the cosmic time t and the radial distance r. At low 
redshift, it can be described by a Lemaitre-Tolman-Bondi 
(LTB) spacetime 0] that depends on two arbitrary func- 
tions of r. Interestingly it was shown that a LTB universe 
reproducing the luminosity distance-redshift relation ob- 
served can be reconstructed without introducing any new 
form of matter [l6| and it was shown that homoge neity 
cannot be proven using only background quantities [17J. 
This implies that the Copernician pri ncip le has to be 
tested observationally (see e.g. Refs. (T^j) as much as 
possible, particularly under the light of some new pro- 
posals [l9| . 

We will not consider this interesting possibility in the 
following and we assume that the universe is well de- 
scribed by a Friedmann-Lemaitre spacetime. From a cos- 
mological point of view, the dynamics of the background 
expansion is characterized by the matter content of the 
universe, that is a list of fluids 2 with their equation of 
state. All background observations are then related to 
the function H(o)/Hq where Ho is the Hubble constant 
today. The extra degrees of freedom, often referred to 
as dark energy and needed to explain the data, can be 
introduced as a new kind of matter or as a new property 
of gravity. 

Let us clarify this. General relativity relies on two 
assumptions: (1) gravitational interaction is described 
by a massless spin-2 field and (2) matter is minimally 
coupled to the metric which, in particular, implies the 
validity of the weak equivalence principle. It follows that 



S = 



16nG 



-gd 4 x + S m [ma.t;g tlI/ ] 



(1) 



2 As long as the background dynamics is concerned, the symmetry 
imposed by the cosmological principle is such that only perfect 
fluids can be considered. This is why the parameterization of 
dark energy reduces to an equation of state at this level. Con- 
sider a comoving observer with 4- velocity « M perpendicular to the 
hypersurfaces of homogeneity. The line element of the FL space- 
time takes the form ds 2 = — (n M dx M ) 2 + r y IJ/l/ dx )1 dx 1 ' with 7 M „ = 
+ ti(i«u. As a symmetric rank-2 tensor, the stress-energy 
tensor of any matter should decompose as = Au^Uy + B^^ V 
and A = T^u^u" reduces to the energy density measured by 
this comoving observer and B to the isotropic pressure. Thus, 
there is no extra-assumption at this stage. 
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where g^ is the metric tensor and R its Ricci scalar. 
S m [mat; g^] is the action of the matter fields. Each class 
of models will account for a modification of the minimal 
action and our discussion restricts itself to models rely- 
ing on a field theory. This has the advantage to discuss 
which are the new degrees of freedom and to which ex- 
tent the theory is well defined. It also implies that we do 
not consider models in which the Friedmann equations 
are modified without an underlying field theory which 
we consider on the same footing as an ad hoc parameter- 
ization of the equation of state. 

In the first approach one assumes that gravitation is 
described by general relativity while introducing new 
forms of gravitating components beyond the standard 
model of particle physics to explain the observed accel- 
eration of the universe. This means that one adds a 
new term Sd e [ip; g^v] in the action while keeping the 
Einstein-Hilbcrt action and the coupling of all the fields 
(standard matter and dark matter) unchanged. 

The other route is to allow for a modification of grav- 
ity. This means that the only long range force that 
cannot be screened is assumed not to be described by 
general relativity. Various ways to extend the minimal 
action JTJ have been considered, modifying the Einstein- 
Hilbert action or the coupling of matter. Whatever the 
modification, one has to introduce new degrees of free- 
dom (not necessarily scalar) and even though we refer 
to these models as a modification of gravity, we have to 
keep in mind that they involve also new matter. 

This being defined, we can distinguish four general 
classes ^jf 3 . First are models in which gravity is not 
modified 

1 . Class A consists of models in which the acceleration 
is driven by the gravitational effect of new fields 
which are minimally coupled to gravity. It follows 
that they are not coupled to the standard matter 
fields or to dark matter and that one is adding a 
new sector 
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FIG. 1: Summary of the proposed different classes of mod- 
els. As discussed in the text, various tests can be designed to 
distinguish between them. The classes differ according to the 
kind of new fields and to the way they couple to the metric 
(? MJ/ and to the standard matter fields. Left column accounts 
for models where gravitation is described by general relativity 
while right column models describe a modification of gravity. 
In the upper line classes, the new fields dominate the mat- 
ter content of the universe at low redshift. Upper-left models 
(class A) consist of models in which a new kind of gravitating 
matter is introduced. In the upper-right models (class C), a 
light field induces a long-range force so that gravity is not de- 
scribed by a massless spin-2 graviton only. This is the case of 
scalar-tensor theories of gravity. In this class, Einstein equa- 
tions are modified and there may be a variation of the fun- 
damental constants. The lower-right models (class D) corre- 
sponds to models in which there may exist massive gravitons, 
such as in some class of braneworld scenarios. These models 
predict a modification of the Poisson equation on large scales. 
In the last models (lower-left, class B), the distance duality 
relation may be violated. 



models 23] which try to unify dark matter and dark 
energy, and K-essence models |2 ll models invoking 
scalar fields with a non-canonical kinetic term. 



Sdelde;^] 

to the action . To explain the acceleration of the 
universe, the new matter must have an equation of 
state, w = P/p smaller than —1/3. As standard 
examples of this class of models, we can cite the 
standard A CDM model (which has the property to 
involve only 1 new parameter and no new field), 
quintessence models pfl| which invoke a canoni- 
cal scalar field slow-rolling today, solid dark matter 
models [2lJ induced by frustrated topological de- 
fects networks, tachyon models [22], Chaplygin gaz 



3 A similar classification can be applied to dark matter issue. The 
list of models cited here is indeed incomplete and we are not 
aiming at an exhaustive review of specific models. See Ref. 0, 
130 for that purpose. 



2. Class B introduces new fields which do not domi- 
nate the matter content so that they do not change 
the expansion rate of the universe. These fields are 
coupled to photons and thus affect the observed 
dynamics of the universe but not the dynamics it- 
self. In particular they are not required to have an 
equation of state smaller than —1/3. An example 
is provided by photon- axion oscillations in an ex- 
ternal magnetic field |2Jj| which aims at explaining 
the dimming of supernovae, not by an acceleration 
expansion but by the fact that part of the photons 
has oscillated into invisible axions. In that par- 
ticular case, the electromagnetic sector is modified 
according to 

A specific signature of these models would be a vi- 
olation of the distance duality relation [l3l l2rj| , a 
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possible violation of the variation of CMB temper- 
ature with redshift [T~ oc (1 + z)] which seems to 
hold observationnally [2]} and in the future, the 
determination of the luminosity from gravitational 
waves will be insensitive to such a coupling. 

Then come models with a modification of gravity. 
Once such a possibility is considered, many new mod- 
els arise |2^. Indeed in considering these modifications, 
one needs to take great care about the fact that the new 
theory is well defined [2!| and stable both outside and 
inside matter (as an example theories such as /(R^) or 
f{Rfj.v P a) involve extra massive spin-2 ghosts [20]; we also 
stress the Ostrogradksi theorem [3lJ as recently exposed 
in Ref. |2^). We can distinguish two classes, 

3. Class C includes models in which a finite number of 
new fields are introduced. These fields couple to the 
standard model fields and some of them dominate 
the matter content (at least at late time). This is 
the case in particular of scalar-tensor theories in 
which a scalar field couples universally and leads 
to the class of extended quintessence models [32j . 
chameleon models j3j| or f(R) models depending 
on the choice of the coupling function and potential. 
For these models, one would have a new sector 

and the couplings of the matter fields will be mod- 
ified according to 

S m [mat;g^ v ] -> 5 m [matj; A?(tp)g MV ]. 

In the case where the coupling is not universal, a 
signature may be the variation of fundamental con- 
stants [34| and a violation of the universality of free 
fall. This was argued to be a general signature of 
quintessence models [35| . This is the case in partic- 
ular in the runaway dilaton model [3(|. This class 
also offers the possibility to enjoy w < — 1 |37l l3S| 
with a well-defined field theory and includes models 
in which a scalar field couples differently to stan- 
dard matter field and dark matter [39l |. 

4. Class D includes more drastic modifications of 
gravity with e.g. the possibility to have more types 
of gravitons (massive or not and most probably 
an infinite number of them). This is for instance 
the case of models involving extra-dimensions such 
as e.g. multi-brane models |4£j, multigravity EL 
brane induced gravity or simulated gravity jig • 
In these cases, the new fields modified the gravita- 
tional interaction on large scale but do not neces- 
sarily dominate the matter content of the universe. 
Some of these models may also offer the possibility 
to mimick an equation of state w < — 1. 

These various modifications can indeed be combined to 
get more exotic models. Since models of classes C and 



D involve departure from general relativity, they have to 
pass the local tests on the deviation from general relativ- 
ity (see Refs. [37l EH for examples) both in weak field in 
the Solar system and in strong field (timing of pulsars). 
Both classes lead to a modification of the Poisson equa- 
tion on sub-Hubble scales that can be tested using weak 
lensing or the large scale structure [45| . 

Let us stress two important points. First, these models 
do not address the cosmological constant problem (see 
Refs. 0[I3 for reviews on this aspect). Second, at the 
moment the analysis of the various cosmological data do 
not push for any time dependent equation of state and are 
completely compatible with the standard ACDM model 
(see e.g. Ref. Q for a general review on the constraints, 
Refs. PI 0,|5( for the constraints on a constant equation of 
state 4 and Ref. for the case of quintessence models) 
and there is no need at present for anything else but a 
cosmological constant [43 . 

In order to have some handle on the physics behind 
dark energy, one needs to detail the various signatures 
of these models. They may be of different natures: (1) 
dynamics of the background (related to the equivalent 
equation of state for each model), (2) properties of the 
growth of structure on sub-Hubble scales, (3) properties 
of perturbations involving super-Hubble scales, such as 
the CMB, (4) non-linear clustering, (5) local tests and 
(6) strong field effects such as gravitational waves pro- 
duction and low acceleration regimes (e.g. galaxy rota- 
tion curves ). 

From an observational point of view, we can briefly 



4 Let us make a short parenthesis concerning models with a con- 
stant equation of state. Such models are often used to evaluate 
the deviation from a ACDM (see e.g. Refs. QQS)- Can we 
determine what kind of models we test with this assumption? In- 
deed, using the standard perturbation equations means that we 
are dealing with class A models. We can use Eqs. I39I4UI of sec- 
tion III.B.2 to characterize the quintessence models sharing this 
property. If w^ c is constant then pde = Pdeo(l + ,z) 3 ( 1+1 " do ' so 
that, restricting to a fiat universe, E 2 = O m o(l + z) 3 + f2d e o(l + 
£)3(l+m d e). With the same notations as above, one gets 



and 



V(z) 



3w de 



■ ™dc 



fWl + z) 3(1+1 " do) 



arg sinh 



- arg sinh 



(1 + Z ) 3w de 



^dcO 



Clearly using a constant equation of state (and larger than —1) 
without modifying the perturbation equations tests only a very 
specific class of potential and its meaning in terms of physical 
models is far from being clear when one allows w dc < — 1. 
Even though we do not aim at designing models explaining both 
the acceleration of the universe and galaxy rotation curves, the 
effect of the new fields and in particular the modification of grav- 
ity, if any, has to be taken into account in that regime and may 
have some implication for dark matter. 
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summarize the data sets that allow to probe these various 
regimes. (1) The background dynamics can be probed 
from the Hubble diagram [and more particularly using 
SN la; the angular distance-redshift relation from strong 
lcnsing; comoving volume; number counts vs. redshift 
using galaxy number counts; baryon acoustic oscillation 
(BAO) 6 ; the CMB shift parameter; big-bang nucleosyn- 
thesis; and in the future the determination of luminosity 
distance from gravitational waves]. (2) The growth of 
structures in the Newtonian regime can be obtained from 
3D-galaxy surveys (luminous matter), halo number count 
vs. redshift using clusters of galaxies (X-ray, optical/NIR 
counts and velocity dispersion, Sunyaev-Zel'dovich effect, 
strong and weak lensing) and galaxy redshift distribution 
(dark halos), Lyman-a forest, weak lensing (2D- survey 
but with the possible future development of tomography) 
that gives access to the matter and gravitational poten- 
tial power spectra at different redshifts and on different 
scales. It is important to separate the linear and non- 
linear regimes that involves the regime (4). (3) Super- 
Hubble perturbations affect mainly the CMB in general 
and the integrated Sachs- Wolfe which, in particular when 
correlated with large scale structures, contains some in- 
formation on dark energy. (5) Local tests are usually 
performed using the parameterized post-Newtonian for- 
malism in the Solar system and (6) concerns mainly the 
timing of pulsars and the emission of gravitational waves. 

From a theoretical point of view, it is important to de- 
termine which models are in fact two versions of the same 
model which means to determine the nature and cou- 
plings of the degrees of freedom (e.g. the fact that f (R) 
models are nothing else but scalar-tensor theories [43 
or that theories involving f(R, OR, . . . , □™i?) are (n + 1) 
scalar-tensor theories |50| ; the relation between k-essence 
and quintessence |U) and to which extent the same set 
of data can be reproduced by various theories. In this 
sense the determination of the number and nature of the 
new degrees of freedom are important and gives a the- 
oretical estimate of the complexity of the theory much 
more accurate than the number of free parameters of a 
model. Obviously, a subset of all available data can be 
explained by models in different classes (see below) so 
that our classification will also help quantify the habiliy 
of a given set of data to distinguish between theoretical 
models, which is different from determining which family 



The wavelength of the BAO is mainly determined by the 
wavevector k = 2tt/s with 




c a (z)[H(z)/H Q ]- 1 dz 



where c s = 1/ J 3(1 + 3pb/4p7)- Even though it relies on the dy- 
namics of the photon-baryon fluid on sub-Hubble scales, it would 
be interesting to check that modification of gravity does not alter 
this picture. Let us also emphasize the existence of distortions 
due to the spacetime shear which enable a possible test of the 
Copernician principle using the Alcock-Paczynski test l48l . 



fits best. 

From a more phenomenological point of view, it is of- 
ten usefull to rely on an effective parameterization of the 
equation of state to describe the change in the dynamics 
of the expansion. Although it is a key empirical informa- 
tion on the rough nature of dark energy, a detailed de- 
scription of its properties demands more thoughtful data 
interpretation. To be useful, the parameterization has 
to be realistic, in the sense that it should reproduce pre- 
dictions of a large class of models, it has to minimize 
the number of free parameters, and it has to be related 
to the underlying physics (see e.g. Ref. [12] )• Because 
the result of a data analysis will necessarily have some 
amount of parameterization dependence 53], choosing 
a specified physical model strategy seems preferable to 
break degeneracies. In particular, it enables to compute 
without any ambiguity their signature both in low and 
high redshift surveys, such as the cosmic microwave back- 
ground. The increasingly flourishing number of models 
hampers to provide a comprehensive set of unambiguous 
predictions to constrain physical models one by one with 
present-day observations, but there are still several bene- 
fits in exploring dark energy this way, in particular when 
weak lensing surveys are used together with CMB obser- 
vations. In this respect also, our classification may turn 
to be usefull. 

It follows that one is stuck between pragmatism of the 
data analysis favouring parameterization of the equation 
of state (bottom-up approach) and the theorist point of 
view (top-down approach). It is indeed clear that the 
physical nature of dark energy will not directly come out 
of the observations and we need obviously to go beyond 
the measurement of w to have any idea of the models 
or classes of models that are likely to be good candidate 
both to explain the observed universe and from a theo- 
retical perspective. In this respect, constructing target 
models in each classes is a key issue. Both approaches 
are indeed complementary. 

Here, we want to illustrate how perturbations are of 
importance when one wants to tackle down the physics 
behind the acceleration of the universe. In Section [HI 
we start by recalling the properties of the background 
evolution and then we revisit the description of the per- 
turbations and in particular we focus on their properties 
that can shed some light on the class of theoretical mod- 
els that are best suited (see Fig.Q. This will lead us to 
propose a chain of tests to investigate dark energy. In 
Section IlIII we consider a first model inspired by brane 
world cosmology and then turn, in Section IIVI to the 
more involved case of scalar-tensor theories of gravity. 



II. DISTINGUISHING MODELS 

In this section, we compare the information that can 
be extracted from the background dynamics and the dy- 
namics of density perturbation at small redshift. This 
will lead us to propose a chain of tests to characterize 
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the properties of the dark energy. 



A. Background evolution 

Assuming the validity of the Copernician principle, the 
equation of state of the dark energy is obtained from 
the expansion history, assuming the standard Friedmann 
equation. It is th us g iven by the general expression (see 
e.g. Refs. 00 E3) 



3fla c w = — 1 + Q.k + 2g, 
q being the deceleration parameter, 



a a 

"^2 



„ 1 N dlnff 2 

-1 + ^1 + *—, • 

2 dz 



(2) 



(3) 



This expression does not assume the validity of general 
relativity or any theory of gravity but gives the relation 
between the dynamics of the expansion and the property 
of the matter that would lead to this acceleration if gen- 
eral relativity described gravity. Thus, it reduces to the 
ratio of the pressure to energy density only under this 
assumption. All the background information about dark 
energy is encapsulated in the single function w. 

As long as the background dynamics is concerned, all 
the observations are related to the function E{z) defined 

by 



In particular, the angular distance is given by 



D A (z) = D Ho - 



1 



(4) 



(5) 



where Du = c/Hq. The radial distance \ is defined by 

1 [ z dz' 



aoH J E(z') 
and the angular diameter distance Sk is given by 

Sk(x) 



sm(VK X )/VK 
X 

sinh(V— Kx)/ 




(6) 



(7) 



In the case where the distance duality relation holds, the 
luminosity distance is given by 

D L {z) = {l + zfD A (z), (8) 

but this may not be the case in models of the class B. 

B. Parameterization of the equation of state 

Most data, and in particular supernovae data, are 
being analyzed using a general parameterization of the 



equation of state. These parameterizations are useful to 
extract model-independent information from the obser- 
vations but the interpretation of these parameters and 
their relation with the physical models are not always 
straightforward. 

As a first example, recall that general parameteriza- 
tions of the equation of state as 

w(a) = w(a ) + [w(a m ) - w(a )} T(a, a t , A) (9) 

were shown to allow an adequate treatment of a large 
class of quintessence models [5^, |5?j. It involves four 
parameters {w(ao), w(a m ), at, A} and a free function T 
varying smoothly between one at high redshift to zero to- 
day. Even though it reproduces the equation of state of 
most quintessence models, it is not economical in terms 
of number of parameters since most quintessence poten- 
tials involve one or two free parameters. If one assumes 
that the parameterization is supposed to describe the dy- 
namics of a minimally coupled scalar field, the knowledge 
of w is indeed sufficient but in a more general case one 
would need more information. The background dynamics 
depends on the potential and its first derivative, which 
can be related to w and its derivative. Accounting for 
perturbations, one needs to know the second derivative 
of the potential which can be inferred from id [Hsj . 

Since we expect dark energy to have observable conse- 
quences on the dynamics only at late time, one can con- 
sider an equation of state obtained as a Taylor expansion 
around a pivot point, 



1 



1 + z* 



1 



(10) 



This form depends on only three parameters and is a 
generalization of the parameterization proposed by [59| 
and then [60j where z* = 0. Whatever the expression for 
w, it is easy to sort out that 



Pdc = -3(1 + w)p dc 



(11) 



where a prime denotes a derivative with respect to the 
number of e-folds, p = ln(a/ao). It follows that 



pde = Pdco exp 



(1 + w)dp 



and thus 
Hi 



£! m oe 



-3p 



H, K0 e 



-2p 



£!dco e 



-3/ p (l+S)dp 



(12) 



(13) 



Two considerations are in order when using such a pa- 
rameterization. First, the redshift band on which this is 
a good approximation of the equation of state is a priori 
unknown. Clearly, compared with the form (O, it is un- 
likely to describe dark energy up to recombination time. 
Secondly, when combining observables at different red- 
shift such as weak lensing, Sn la and CMB, one should 
choose the value of z* in such a way that the errors in w* 
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and w a are uncorrelated 61J and the pivot redshift is the 
redshift at which w is best constrained. In particular, it 
was argued that it is important to choose z* for distance- 
based measurements. The problem lies in the fact that 
the pivot redshift is specific to the observable. In this 
respect, dark energy models defined by a Lagrangian are 
more suitable, yielding to a definite equation of state as 
a function of redshift, hence more general than a Tay- 
lor expansion around a pivot point. Eventually, one can 
read out the values of w* and w a at whatever redshift. 
This complication, arising when one wants to combine 
datasets with different z», will also make it more diffi- 
cult to infer constraints on the physical models from the 
constraints on the parameterization (see Refs. for 
examples) . 

There is an alternative way to get a first hint on the na- 
ture of dark energy. It may be useful to consider the plane 
(w,w ! ) where w' = dw/dp. It was recently shown |63tl64| 
that quintessence models occupy a narrow part of this 
plane (assuming z* — 0). As discussed later on, the way 
to relate these parameters to a physical model (that is 
the reciproque) is difficult. 

In conclusion, a class of equations of states or the 
(w,w') analysis are usefull tools to constraint a class of 
theoretical models independently of the details of each 
model and the parameterization must be designed for 
this class of models. In such a case the position in the 
(w, w') can also be related to a particular model in this 
class. For instance, in the case of quintessence models 
where the acceleration is due to the gravitational effect 
of a slow-rolling minimally coupled scalar field, the posi- 
tion in the (w,w') plane can be related to the slow-roll 
parameters 

characterizing the shape of the potential V close to the 
value of the scalar field today since w ~ — 1 + 2e/3 and 
w' ~ — 4e(?7 — 2e)/3. Unfortunately, this cannot deter- 
mine the model since the same position can be reached 
by various models from different classes. 

For background evolution, we thus have various pa- 
rameterizations that describe class A models and include 
the ACDM as a subcase. They enable to quantify, within 
this particular class, how close from a ACDM the back- 
ground dynamics allows to be. In the particular case of 
quintessence, the relation with the physical parameters 
is clear but one should impose a prior on the parameter- 
izations to ensure that w > —1 as is the case for these 
models (see e.g. Ref. for a comparison of a data anal- 
ysis based on a model and a parameterization). Indeed, 
great care is required when interpreting constraints ob- 
tained from such a parameterization when relaxing this 
prior. 



C. Taking perturbations into account 

From the study of the background dynamics, one can, 
in principle, determine or constrain w(a). Other sets of 
data, such as weak lensing and large scale structures, 
involve density perturbations. Assuming that the metric 
of spacetime takes the form 

ds 2 = -(1 + 2$)dt 2 + (1 - 2V)a 2 j ij dx i dx j , (15) 

where Yy is the metric of the spatial section, the evo- 
lution of density perturbations on sub-Hubble scales is 
dictated by 

S m + 2H6 m - 4^A$ = S dc . (16) 

This equation derives from the Euler and conservation 
equations on sub-Hubble scales, that is from the con- 
servation of the stress-energy tensor of matter. It could 
inherit a source term Sdc =/= if the presureless matter is 
coupled to other matter species (e.g. in class C). 

To close this equation, one needs to use Einstein equa- 
tions to express $. Assuming general relativity and a 
ACDM model, the Poisson equation takes its standard 
form 

A* = AnG Pm a 2 6 m = ^l m H 2 a 2 5 m , (17) 

and the two gravitational potentials are related by 

$ - * = (18) 

since the anisotropic stress of radiation is negligible in 
that regime. When considering dark energy models, one 
has to allow possible modifications of these two equa- 
tions. 

First, the second Einstein equation can enjoy a non- 
vanishing anisotropic stress due to the dark energy sector 
so that it takes the general form 

A($ - *) = IT d0 (19) 

where n<j e can be decomposed as Ildc = 87rGPdoa 2 A7rd . 
This term leads to a general source term Sd e + Hde /a 2 
m Eq. O but it is worthwile making the distinction 
because alone will be involved in lensing observations 
that depend on $ + \&. All classes can induce such a term. 

Second, the density perturbation of the dark energy 
may be non negligible in the Poisson equation and, as in 
the case with a modification of gravity, one can expect a 
modification of proportionality coefficient that may even 
be scale dependent. So we will assume that the general 
form of the Poisson equation is, in Fourier space, 

-fc 2 $ = 4nGa 2 F m (k, H)p m S m + A do . (20) 

If F m depends on k then the modification will not be 
degenerate with the normalisation of the amplitude of 
the matter power spectrum and can be tested [65j other- 
wise one will have to use the growth history to determine 
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Validity of 

Copernician Principle 



No Inhomogeneous 
*~ spacetimes 



lYes 



Friedmann-Lemaitre spacetime 
a(t); E{z) = H/Hq 



Background dynamics 
D L (z); D A (z),... 



D L (z) = (l + z) 2 D A (z)l 



No 



, Class B 



jYes 

SfldcW = —1 + Qk + 2q 



t-,v • , jit Yes Fit background 
Physicajmodel ? ► + pertur S bations 

Design a parameterization -* 

I 

Fit background [CI, w(z), . . .] 
Information from perturbations 



I 



Test of Poisson equation 
D+ = D + ? 



Yes 



|No 

New tests (constants,. . . ) 
modification of gravity 
Classes C and D 



Class A ■ 



w = —I? 

|Yes 

ACDM 



whether F m = 1 or not. Such a term can arise in classes 
C (see e.g. Ref. [66|1 and D (see Ref. We emphasize 
that for classes A and B, the growth factor is scale in- 
dependent and that the primordial scale-dependence and 
the growth factor decouple. This may not be the case in 
classes C and D when F m is scale dependent. Thus, it is 
safer to distinguish P m (k, z) and P$(fc, z). 

It follows that in general, the equation of evolution of 
the growth factor is expected to take the form 

D + 2HD-4irG P F m D=^+S dc . (21) 

Indeed, such an equation is not closed and one needs 
either to propose a parametrization of the source term 
or explicit the evolution equations of the dark sector. 
Note that while the growth factor depends on (-Fm, Ad e , 
6>d c ), weak lensing involves the combination A($ + = 
8irGpa 2 F m 5 m + 2Ad c — Rde [68| so that combining them 
may be very fruitful. 

In conclusion, the evolution of perturbations at low 
redshift of the various classes of models can be charac- 
terized in the Newtonian regime by (£>dc, Ade, -Fm, Bd c )- 
Models of the class A satisfy (iSde, Ad e , F m , Rde) = 
(0, Ade, 1, Ilde); for instance, a ACDM model has 

Ilde) = (0,0,1,0) (22) 

while for quintessence models, <5pdc = ipStp — tp 2 & + V'Sip 
so that 

(S dc , A dc , F m , Ilde) * (0,0,1,0) (23) 

where ip is the quintessence field and Sip its perturbation. 
This shows in particular that quintessence and ACDM 
can only be distinguished, at low redshift, on the basis 
of their equation of state. Models of the class B have 
the same characteristics than a ACDM for the density 
perturbation evolution but have an effective equation of 
state that may differ from w = — l 7 . Models of the classes 
C and D have more general (Sd e , Ade, F m , Tide) and in 
particular F m ^ 1. We shall give two examples further 
on. 



D. Summary 



We have argued that the different models appearing in 
the classification presented in Fig. ^ can be characterized 
by a set of functions including the equation of state to 
describe the background dynamics and (<Sde, i5de, -F m , Bide) 
for the evolution of the perturbations at low redshift. In 



7 Note that in class B, the effective equation of state derived from 
the observed dynamics of the universe is different from the "true" 
equation of state which may be obtained from the growth factor. 



FIG. 2: A possible chain of tests to unveil the nature of dark 
energy. The goal is to start from the more general hypothesis 
and to incorporate new data and information one by one in 
order to check at each step if the hypothesis underlying the 
equations used in the analysis hold or not. Here D+ refers 
to the growth factor predicted from w(z) assuming general 
relativity. In particular, it may turn out that a deviation from 
the Poisson equation may be detected while no deviation from 
w — — 1 is established. This would however require to extend 
the minimal ACDM. 

most cases, the dependence on the primordial spectrum 
and on the dark energy sector decouple in the Newtonian 
regime. Indeed, more information would be needed to 
treat CMB anisotropies and strong field phenomena. 

This description can help to get some handle the na- 
ture of dark energy and in particular can give a way to 
organize the theoretical landscape and the way to distin- 
guish between classes. In particular, it will be usefull to 
define both a target model and parameterizations in each 
class so that the distance to the ACDM in each category 
can be quantified. Such parameterizations exist for the 
class A but they need to be generalized to other classes 
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and the constraints on the parameters of these parame- 
tcrizations for them to describe physical models actually 
in these classes should be established. 

Let us also mention a difference with a simple fit of 
models to data. In this case, one ususaly minimizes the 
information of the data and the theory together, I(D, T) 
taking into account a penalty, I(T), that de pen ds on the 
number of the free parameters of the theory [69j . The to- 
tal information reduces to I(D, T) = I(D\T) + I(T) with 
I(D\T) being the opposite of the likelyhood and various 
prescriptions for I(T), all depending on the number of 
free parameters, exist but in fact it depends on the com- 
plexity of the theory. It is in general difficult to define 
but it is related to the number, nature and couplings of 
the new fields of the theory. To this respect, the ACDM 
model is the simplest since it involves only 1 new param- 
eter and no new field. Relying on null test (or exclusion 
tests) is a way to indicate the degree of complexity re- 
quired without having to consider such a statistical issue 
in the first place (see Ref. [7(j for the implication con- 
cerning dark energy). 

It follows from the previous discussion that various 
questions arise in order to grasp some of the physics be- 
hind dark energy. We propose the chart in Fig. [21 to 
address them and outline the dark energy sector. It is 
based on a series of consistency checks trying to exclude 
the most simple models in order to determine how com- 
plex a minimal model should be. It also tries to determine 
to which extent data can discriminate between the vari- 
ous classes and calls for the construction of new tests. At 
the heart of it, is the question of the necessity to extend 
the standard ACDM and to determine, in a robust way, 
how close from it we must be. Indeed demonstrating the 
time variabiliy of w is difficult to proove but it may be 
that other signatures will be easier to obtain 71] . 



A. Summary of the cosmological properties 



III. FIRST EXAMPLE: DGP COSMOLOGY 



The Dvali-Gabadadze-Porrati (DGP) brane-world 
model 0j E2 i s constructed from a brane embedded 
in a 5-dimensional Minkowski bulk. It was shown that 
because of gravity leakage to the bulk, it leads to a mod- 
ification of gravity on large scales and could explain the 
recent acceleration phase. The characteristic scale, r c , 
at which the modification of gravity manifests itself is 
related to the two mass scales of the model, the usual 4- 
dimensional Planck mass and the 5-dimensional Planck 
mass, M 5 . When M 5 - 10 - 100 MeV, r c is typically of 
the order of c/Hq. Interestingly, this model involves only 
I extra-parameter, as the standard ACDM, but involves 
new fields. 



In the DGP model (see Refs. |42|, |72j for a description 
of the model), the Friedmann equation takes the form 




TZi+tt- (24) 



and we consider the e = 1 case. r c is the new parameter 
of the model. We focuse to the low redshift universe 
where p is dominated by presureless matter. We define 

I 



^cO 



so that it rewrites as 



e 2 + ^o(i + ^) 2 



where we have defined 

&(z) = y/n^ + ^n c0 + n m0 (i + zf. 

It follows that we have the constraint 

e 2 (o) + n KO = (v^ + v /n co + fimo) 2 + 



the solution of which is 

(-l + fi m0 



n. 



c0 



4(1 -SIko) 



(25) 

(26) 

(27) 

Qko = 1 
(28) 

(29) 



By comparing Eq. (|2(j[) with its standard form, we de- 
duce that 



87rG/9dc 

or equivalcntly 

Od c — 2\/ O 

where 



(30) 



(31) 



^cO 



HE 
H 2 



On 



fi ro0 (l + z) 3 §|. 



In particular, Q^c + + = 1 is equivalent to the 
constraint (|28|l when evaluated today. Using the expres- 
sion H26|) we obtain that 



1 + w(z) = 



fimo(l + zf 



(32) 



(33) 



2e(z) v /fi c o + fi m o(l + z) 3 
whatever the curvature. It is easy to sort out that 

, _ dw _ 3 \/n^Q, m0 (l + z) 3 

W ~ d^ ~ ~4[ft c0 + a„o(l + ^) 3 ] 3 / 2 ' 

As a first consequence of this analysis, Fig. shows 
that today the background evolution of DGP models can 
be mimicked by quintessence models. Note that the zone 
of the plane (w , dw /dp) that corresponds to DGP mod- 
els remains very thin even when we allow Qko to vary 
between —1 and 1 and f2 m o between and 1. 



FIG. 3: In the plane (w,w'), the filled zone corresponds to 
quintessence models while the very thin black zone corre- 
sponds to DGP models. We have allowed Qko to vary be- 
tween — 1 and 1 and fi m o between and 1. 



B. Models with equivalent background evolution 

1. Fitting the equation of state with a parameterization 

The matching to a given equation of state can be per- 
formed in various ways. First, once z* has been chosen, 
we can require to fit (w,w') at the pivot redshift, so that 



FIG. 4: Comparison of the fits of the DGP equation. The solid 
line corresponds to the fit defined in Eq. I|34jl while the dashed 
and dotted lines correspond to the fit defined in Eq. 13511 re- 
spectively with z p = 3 and z v — 1000. We have assume 
Omo = . 3 and flap — . 




1/(1+ z) 



-(1 + 2.XW- ( 34 ) 



FIG. 5: Comparaison of linear growth factor for the same fits 
as in Fig.EJto the equivalent quintessence model. 



Another route is to require that the equation of state has 
its correct value today and at a redshift z p so that 



w(z*), w a = [w(z p ) - w(z*)] 



(1+Zp)(l + Z«) 

(35) 

In Fig.QJ we compare these two fits of the DGP equa- 
tion of state (|32|1 assuming z* — and choosing z p = 3 
and z p — 100. By construction both matching lead to the 
correct value of the equation of state today (w* = Wo) 
but the derivative differs by 6.6% and —5.3%. Indeed 
these are very small deviations but they may be impor- 
tant while trying to infer whether a constraint on (wo, w' ) 
points toward a DGP model or a quintessence model. So 
the question that arises is how a measured (u>*, w a ) is ac- 
tually related to the constraint plane depicted in Fig. [3] 
and to the physical parameters of the model. 

Interestingly, as shown in Fig. \5\ the error induced 
on the growth factor coming from the error in the fit- 
ting of the equation of state is smaller than 1%. This 
implies that one cannot distinguish between the three 
parametrizations studied in Fig.^Jfrom the study of per- 
turbations. 



2. Equivalent scalar field model 

We can even go a bit further and construct explicitely 
the quintessence model that reproduces the DGP dynam- 
ics for the background. 

Quintessence models reduce to the dynamics of a min- 
imally coupled scalar field, Q, evolving in a potential 
V(Q). The model is completely specified once the poten- 
tial is given. It is thus clear that the energy density is 
given by 



1 .. 
PQ = t;Q 



V 



and the equation of state is 



wq 



Q 2 -2V 
Q 2 + 2V' 



(36) 



(37) 



Inverting these two relations, one can reconstruct the po- 
tential [Z3 mimicking a dark energy component 8 charac- 



See also Refs. 180.1 for a similar exercice. 
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. 3 

Q-Qo 



FIG. 6: Reconstruction of the potential of a quintessence 
model that gives the same background evolution than the 
DGP model when Q m o = 0.3 and Qko ~ (solid), Qko ~ 0.3 
(dashed) and £Iko ~ —0.3 (dotted). 



terized by {p de (a),w dc (a)} since we have 



Q 2 = Pdc(a)[l + Wdc(a)], 
V = -Pde(a)[l - Wde(a)]> 



(38) 



It follows that the parametric expression of the potential 
is 



V{ z ) = ^Pde(z)[l ~ W de (z)], 



dQ _ y/p do (z)[l + w dc (z) 
dz 



(l + z)H{z) 



(39) 
(40) 



We can arbitrarily choose the sign of dQ/dz. Us- 
ing Eq. (EOJ and defining V = 8nGV/3H$ and Q = 
^/8irG/3Q, this system take the form 



V(z) = N /O c0 [l - w de (z)]G(z), 
dQ _ {2y/n^[l + Wdc (z)]9(z)} 1/2 



dz 



(l + z)E(z) 



(41) 
(42) 



Once (f2 m o, Qko) is fixed, fl c o is determined and the DGP 
model we want to mimick is completely specified. It fol- 
lows that u> do is given by Eq. and E(z) by Eq. 
and we have a unique potential obtained from the inte- 
gration of the system 1|4H42I1 . 

Figure depicts the potentials of the quintessence 
models that give the same background dynamics that 
DGP model for fi m o = 0.3 and Qko = 0,±0.3. Indeed, 
to choose one or the other model (DGP or quintessence) 
has some importance and takes some theoretical preju- 
dices into account. 

Note also that, as long as we assume that the growth 
of density perturbations, is not modified and is given by 
a form l|23|l . the two models cannot be distinguished. 



C. Growth of density perturbation 

The previous result assumed that the growth of density 
perturbations is dictated by the standard equation. For 



the DGP model, this is not the case and two routes have 
been followed. In the first, five dimensional effects were 
neglected [74], which we will refer to as DGP-4D, so that 



{S de , Ade, F m , lido) — 0,0, 



2Hr n 



.0 



2Hr r - 1 



Unfortunately this setting is not compatible with the 
Bianchi identity. It was recently argued [§3 that when 
five dimensional effects are taken into account we should 
have 

(S dc , A dc ,F m ,II dc ) = ^0,0, 1 + ^^p m a 2 5 D 
with 




(43) 



We shall refer to this case as DGP-5D. The influence of 
this modification of the growth factor on lensing observ- 
ables was recently studied in Ref. |75| . 

From these definitions, we deduce that in the DGP-4D 
case, the equation of evolution of the growth factor 121|) 
can be rewritten in terms of the growth variable, g = 
D/a, as 



g" + [4+(ln£)V 



+ 



3+(ln£)'--n m — 



E 



'QcO 



5 = (44) 



where we have used that 2Hr c = Ej V^co- In the DGP- 
5D case, it takes the form 



g" + [4+(ln£)V 



3 + (ln£)'--fi E 



1 



1 

3^ 



.9 = 0. (45) 



These two solutions have to be compared to the equiva- 
lent quintessence model with same background dynamics 
for which 

(«S de ,A de ,F m ,II de ) ~ (0,0,1,0). 

Fig. [7| compares the linear growth factor for the DGP- 
4D, DGP-5D and equivalent quintessence model. First, 
we see that the equivalent quintessence model can be 
discriminated from the DGP model on the basis of the 
perturbations (almost 20% difference at redshift 0). This 
clearly illustrates the (trivial) fact that perturbations en- 
code extra-information than the one given by the bak- 
ground dynamics. 

We also see that the "theoretical" uncertainty on the 
way to deal with perturbation can lead to uncertainty 
of order 5-10% on the growth factor. It is important to 
try estimate this uncertainty, particularly for models that 
deviates significantly from the standard ones. 
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FIG. 7: (Top) Comparison of the linear growth factors of the 
DGP-5D (solid), DGP-4D (dash) and equivalent quintessence 
(dot) models. (Bottom), the relative error with respect to the 
DGP-5D have been computed and can become large. We have 
assumed f2 m o = 0.3 and £Iko = 0. 



D. Summary 

We have considered a DGP model that belongs to the 
class D. Using the fact that the effective equation of state 
is explicitely known we have compared it to a standard 
parameterization. We have constructed a quintessence 
model (class A) sharing exactly the same background dy- 
namics. After sorting out the form of (Sde, Ad e , F m , lido), 
we have shown that knowledge of the density perturba- 
tion can, as expected, help distinguishing these models. 

This gives an explicit example in which the background 
dynamics cannot discriminate between a class A and class 
D models. It also shows that one should also try to quan- 
tify the "uncertainty" of the theoretical models, a diffi- 
cult task indeed. 



IV. SECOND EXAMPLE: SCALAR-TENSOR 
THEORIES 



A. Overview 

In scalar-tensor theories of gravity, gravity is mediated 
not only by a spin-2 graviton but also by a spin-0 scalar 
field that couples universally to matter fields (this ensures 
the universality of free fall). In the Jordan frame, the 
action of the theory takes the form 



S = 



d 4 a 



16ttG 

+ S r , 



■V=g[F{<p)R- 
[ip; 9nv] 



21%)] 
(46) 



where G* is the bare gravitational constant from which 
we define ft* = 8vtG». This action involves three arbi- 
trary functions (F, Z and U) but only two are physical 
since there is still the possibility to redefine the scalar 
field. F needs to be positive to ensure that the graviton 
carries positive energy. S m is the action of the matter 
fields that are coupled minimally to the metric g^ v . 

It follows that the Friedmann equations in Jordan 
frame take the form 

3f(h 2 + ^J = 8vrG*p+ -3HF + U{47) 

-2F\H-^j = 8nG*{p + P) + Ztp 2 

+F - HF. (48) 

The Klein-Gordon and conservation equations are given 
by 



Z{Cp + 3Htp) = 3F V \H + 2iP + j 



p + 3H(p + P) = 0, 



(49) 
(50) 

where a subscript tp stands for a derivative with respect 
to the scalar field. 

These equations define an effective gravitational con- 
stant 



G 



eff 



G./F. 



(51) 



This constant, however, does not correspond to the grav- 
itational constant effectively measured in a Cavendish 
experiment, 



G c 




F 2 



2ZF + 3F2 



(52) 



an expression valid when the scalar field is massless 9 . So- 
lar system constraints on the post-Newtonian parameters 
imply that F 2 /F <4x 10~ 5 so that they do not differ 
significantly at low redshift. 



As a second example, we consider scalar-tensor theories 
of gravity, which is the simplest example of a model of 
class C. 



9 As we shall see below, in e.g. Eq. 1591 . this is a good approx- 
imation on sub-Hubble scales for models where the scalar field 
also accounts for the acceleration of the universe. This may not 
be the case in more intricate models such as chameleon l33l. 
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B. Sub-Hubble perturbations 



from which we deduce 



The general cosmologic al p erturbations have been 
studied in various articles |6(| |J3 and we concentrate 
to the sub-Hubble regime. 

The conservation equation of the standard matter are 
similar to the ones in general relativity so that S^c = 0. 
In the Newtonian regime, it can be shown that 



* - * = -ftp, 



(53) 



so that non-minimal coupling induces an extra- 
contribution to the anisotropic stress, 



IL 



F. 



fASip. 



The Poisson equation takes the form 
G 

A$ = 4Tr-^pa 2 6 m + A dc . 



(54) 



(55) 



Using the definition l|52l) of the gravitational constant, 
we deduce that 




where the last equality has been drawn using the Solar 
system constraints and should hold at low redshift (see 
e.g. Refs. [6(| I?]} for explicit examples of the redshift 
variation of this quantitiy) . It can be shown |66| that on 
sub-Hubble scales, Ad e = — (F ip /2F)A5ip so that scalar- 
tensor models are characterized by 



!5,,..A,u.F„ : .n,,.! - [u,-^ASip^,-^-ASip). 



Fp f v 
F ' F """7 



To go further, one needs to determine Sip and thus use 
the Klein-Gordon equation which reduces, in that limit, 
to 

(A-U w a 2 )5p = F ip A{<S>~2^) (57) 



j j Z 




Sip = F V A§. 



(58) 



It follows that lido and Ade are directly proportional to 
$ and the Poisson equation is given by 



-fc 2 $ = 47T^f 

F 



k 2 F 2 



{2ZF + AF 2 )k 2 



p m a S n 
(59) 

When U vtp is much smaller than the wavelength of the 
modes, which a good approximation in most models such 
as extended quintessence models, Eq. (|58H reduces to 

A$ ~ 4TTG cavP a 2 8 m (60) 
where G cav is defined by Eq. 

This analysis shows that scalar-tensor models are char- 
acterized by 



(<5>do ^ de ' ^m, ^de) 



0,0, 



F 2 



GcavO ' F + 2F 2 



A$ 



\ 0,^,0 

i ^cavO 



(61) 



F m is indeed time dependent but fc independent. Two 
functions of z characterize scalar-tensor models, the 
equation of state and F m , they are related to the two free 
functions of the theory (F and U) and influence both the 
background evolution and the linear growth of density 
perturbations. 



C. Class C vs class A 

As a first exercice, we can quantify the effect of F m and 
compare a scalar-tensor theory with a model of class A 
with the same background dynamics. For that purpose, 
we need to reconstuct the scalar-tensor theory, which 
can be achieved by using th background equations in the 
form Fl 



d 2 F 

dz 2 

zl 



dlnE 



dz 



fdtp 
\dz 



l + z_ 
6 dF 



dF 
~dz~ 



6 



6F 



(l + z) 2 l + z dz 
2U 



2 dlnE n K0 



F 



2U 



(l + z) 2 E 2 



E 2 



-fimO-Fb 



i + zdz (1 + 2) 2 {i + z) 2 e 2 



l + z F 

6 2 flroO^O — 6-=2-OkO 



(62) 



(63) 



where U = UH 2 and f2 m0 = 8nG cav oPmo/3Ha ~ 
87rG*/9 m o/3-Fb-ffo- As was initially shown in Ref. [j^i the 
knowledge of H(z) and D(z) allow to reconstruct the two 
free functions that appear in the microscopic Lagrangian 



of the scalar-tensor theory. 

The reconstruction can be performed in Jordan frame 
but it is usefull to shift to Einstein frame where math- 
ematical consistency is easier to check. To do so, one 
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perform the conformal transformation = F{(p)g^ v so 
that the potential, V, coupling A, and spin-0 degree of 
freedom, <p* , are related to the Jordan frame quantity by 



dcp) 



F 
F 

3 
4 



1/2 M, 
dlni^Y 



+ 



2F(<p) ■ 



(64) 
(65) 

(66) 



We usually define a = din A/ dtp* in terms of which 
the gravitational constant (|52[1 takes the form G cav = 
G*a 2 (l + a 2 ). The latter equation can be rewritten as 



P_ 

1 + z ~F 



( dtp* x 


)-m 2 - 


F" 


"ldln-E 






2F ~ 


2 dz 




1 dln-E 








+ 1 + z dz 


E 2 



F fi ro0 
Z) F"^ (67) 



It follows that A(tp*) and can also be reconstructed 

parametrically. This is important because if* is actu- 
ally the true spin-0 degree of freedom of the theory [7(J 
and must carry positive energy for the theory to be well- 
defined. This implies that we must have (d(p*/dz) 2 > 0. 
We recall that the spin-2 and spin-0 degrees of freedom 
are mixed in Jordan frame so that the positivity of en- 
ergy does not imply 10 that Z(dcp/dz) 2 > 0. In Einstein 
frame, we have access to A(ip*) and check that it is well 
defined and to V(<p*). In particular, the sign of d 2 V/dip 2 
tells us about the sign of the square mass of ip* and it 
will indicate the existence of an instability of the model 
if it were negative. 

To estimate the possible magnitude of the effects of this 
modification of gravity, let us consider a toy example in 
which 



F(z) = 



4p tanh ( -p 

2 V ozc 



l 



^p. tanh ( 

2 \ dz 



(68) 



With this ansatz, the reconstruction is straightforward 
since the potential can be obtained analytically from 
Eq. (ffi2"]) and (p*(z) is deduced by the integration of 
Eq. Eqs. and can then be used to deter- 

mine G cav (z) that enters in the equation for the growth 
of density perturbations. Two questions have then to 
be considered. First can such a model be realized by 
a scalar-tensor theory and second, how different is the 
growth factor compared to the one of the class A model 
sharing the same background dynamics. 

We consider two examples. In the first one, we as- 
sume that Aq = ±10% and an equation of state of the 
form (|10f> . which will translate in a mild change of G cav 



10 If fact this is the case as soon as 3(d In F/dz) 2 > 4(d^*/dz) 2 
which can happen in perfectly regular situations and ip can be- 
come imaginary while ip* remains well-defined. 



between low and hight redshift. Such a modification en- 
sures that we are safe with time variation of the gravita- 
tional constant (see Ref. [34[ for a review and [71] for con- 
straints arising from BBN). Fig.|S|gives the reconstructed 
scalar-tensor theory, that is {A((p*), V(</?*)}, it also com- 
pares A 2 (z) to G cav (z) and shows that they do not differ 
significantly. Indeed the potentials and coupling func- 
tions may not be easy to justify from theoretical basis and 
are completely ad hoc at this point. Note also that not 
all equations of state are possible to reconstruct (in most 
case the scalar-tensor theory is pathologic but the con- 
clusion mostly depends on the ansatz for G cav (z) so that 
this example is nothing but illustrative at this point). In 
the second example, we have played the same game but 
we have allowed for a more drastic change in the coupling 
function, Aq — 0.4, and we have considered an equation 
of state of the form (|10fl with w* = — 1.1 at z* = and 
w a = 1.6, which gives an explicit example of model in 
which wq < — 1. Interestingly, it can be reconstructed to 
give a well defined theory (without ghosts) that crosses 
the so-called phantom barrier. To finish, we compare the 
growth factor in these three models with those obtained 
in the cosmological model with same background dynam- 
ics but which assume general relativity. Not surprizingly, 
the effect on the growth factor is typically of the order of 
A G (see Fig. H0J|. 

This toy (re) construction tells us that if there are some 
tensions between background and perturbation data, so 
that one is willing to abandon class A models, then a 
reconstruction will indicate whether class B models (at 
least the simple one) are able to lift these tensions and 
offer a framework to interpret consistently all data. It is 
important to keep in mind that one needs to reconstruct 
two free functions so that one needs to determine both 
the expansion history and the growth of density pertu- 
bation observationnaly Otherwise, one must make 

some hypothesis on the two free functions [77ll79l |. It was 
in particular demonstrated that a scalar-tensor theory 
with U = cannot mimick the background evolution of 
a ACDM [?2j and it is obvious that the scalar-tensor the- 
ory corresponding to the background and growth factor 
of a ACDM will be defined by F = 1 and U — constant. 



D. Class C vs class D 

We have shown that the background dynamics of DGP 
models can be mimicked by a quintessence model, but 
that this was no more true at the perturbation level. It 
is obvious that a scalar-tensor theory that can reproduce 
a given {w(z), /3(z)} will mimick a DGP model, both 
for the background dynamics and perturbation evolution. 
Such a scalar-tensor theory, if it exists, will be character- 
ized by 



E{z) =E DGP {z), 
G cav = 1 + l/3/?(z) 
G cav0 1 + 1/3/3(0) 



guGp(z). 



(69) 
(70) 
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FIG. 8: Reconstruction of the function A{tp*) (top) and V(ip*) 
(middle) of a scalar-tensor theory reproducing an equation of 
state of the form < lot with z, = 0, u>* = —0.9 and w a = 1.2 
and a gravitational constant of the form 16811 with zq = 2, 
Sz G = 1 and Aq — 0.1 (solid) or A G = -0.1 (dash). We have 
assumed Q m o = 0.3 and £Iko ~ 0. The lower plot compares 
G C av(^) (black) to A 2 (z) (light). 



FIG. 9: Reconstruction of the function A{tp*) (top) and V(ip*) 
(middle) of a scalar-tensor theory reproducing an equation of 
state of the form < HH with z, = 0, w* = —1.1 and w a = 1.6 
and a gravitational constant of the form 16811 with zq = 1, 
Szq = 1 and Ac = 0.4. We have assumed Q m o = 0.3 and 
^ko = 0. The lower plot compares G CSlV (z) (black) to A 2 (z) 
(light). 



In that case, the reconstruction procedure is the follow- 
ing. First, we can eliminate U from Eq. i|tj3|) to get 



z{^ 

az 



d 2 F 
dz 2 
1 dln£; 
1 + z dz 



dlnE 

dz 



E 2 



1 + z 



F-3- 



dF 



E 2 



-Q m0 F . (71) 



This equation, together with Eq. I|52l) . yields a non-linear 
second order differential equation for F(z) with a source 
term given by Eq. {7UJ. 

Fig. ^2 depicts the redshift evolution of the gravita- 
tional constant of the equivalent scalar-tensor theory (if 
it exists). Let us first stress that the time evolution of 
G C av at z = varies a lot with the parameters of the 



DGP model. The constraint 



dlnG r 



dt 



< 6 x 1CT 12 



(72) 



implies that most DGP models, if interpreted as a scalar- 
tensor theory, are not compatible with Solar system 
experiments, since they should fulfill |d In G cav /d2:|o < 
5.86 x 10 _2 /i _1 . Fig. 1111 summarizes the range of cos- 
mological parameters of the DGP model for which the 
reconstructed scalar-tensor theory, if it exists, is compat- 
ible with Solar system tests on the time variation of the 
Newton constant. 

The reconstruction is more easily performed by using 
the Brans-Dicke representation in which F = tp and Z = 
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0.6 0.8 



FIG. 11: (top) Gcav(z) needed for a scalar-tensor theory to 
mimick a DGP model when fl m o = 0.3 and CIko = (solid), 
Qko = 0.3 (dashed) and £Iko = —0.3 (dotted), (bottom): All 
DGP models which parameters are lying in the white zone do 
not fulfill the constraint l|72|l that will appear on the equiva- 
lent scalar-tensor theory. 



FIG. 10: Comparaison of the growth factor of the scalar- 
tensor models with the models sharing same background dy- 
namics but assuming general relativity, (top) Model depicted 
in Fig. (long dash: A G = -10%, short dash A G = +10%, 
light curve: general relativity); (middle) Model depicted in 
Fig. |5|and (bottom) relative deviation. 



u)(if)/if so that 



G» 2lo + 4 
ip 2u + 3 



(73) 



This representation is well-behaved in the sense that if' 2 
remains positive and the energy condition reduces tow > 
—3/2. Using Eq. (f?3")l to express w as 



2lo 



4 - 3A ffDGp(z)l 

AoffDGpOV - 1 



(74) 



with A = (2wo + 4)/(2cj + 3), Eq. (J7TJ reduces to closed 
non-linear second order equation for if. To performed 
the reconstruction, we need to determine fo, f' and ujq. 
While it is always possible to set fo = 1, we need to 
determine the two others. This can be done by using the 



expressions |28|, [76J, 177J of the PPN parameters 

1 



7 



PPN 



1 



W 



(3 



PPN 



1 

1 = -- 



1 



2) 2 <f/ 



4(2w + 3)(w 
and the expression of the time derivative of G, 



r I \ ^ppn 
u ™ / o \ 



1 



r 



PPN 



(75) 
(76) 
to get 

(77) 
(78) 



We see that the question of the reconstruction of the 
DGP cosmological dynamics by a scalar-tensor theory 
depends on the prediction of the DGP theory in the Solar 
system, an issue still not settled 11 |8lj . 



11 In Eq. 1271 and following describing the background dynamics 
of DGP, flmo has been defined as O m o = 8irGp m o /3Hq without 
questioning the fact that G is the Newton constant, the numeri- 
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TABLE I: Summary of the properties of the models considered here. We emphasize the predictions that one computes in these 
models (Y: yes, N: no, controversial: there is no current agreement at the moment). 



model 


class 


background 


Newtonian 


cosmological 


non-linear 


Solar 


strong 


w < -1 








perturbations 


perturbations 


regime 


system 


field 


possible 


Quintessence 


A 


Y 


Y 


Y 


almost 


Y 


Y 


N 


Scalar-tensor 


C 


Y 


Y 


Y 


N 


Y 


Y 


Y 


DGP 


D 


Y 


Y 


controversial 


N 


controversial 


Y 


Y 



TABLE II: Summary of the comparisons performed in our exercices and of the possibility of two models to share the same 
predictions. (bgd=background; Newt. pert.= density perturbation in the Newtonian regime). 





bgd 


bgd + Newt. pert. 


bgd + Newt. pert. + Solar syst. 


DGP vs quintessence 


Y 


N 


N 


DGP vs scalar-tensor 


Y 


? 


N 



Note that in Ref. |67|. it was suggested that Eq. H45|) 
implies that uj(z) — 3[/3(z) — l]/2. Comparing with 
Eq. I|73|) . this would mean that G cav (0) is identified with 
G*/Fo, which is the case only if ujo ^S> 1 and that the 
redshift dependence of if is neglected. In such a case, 
one can check that lo < —3/2 so that there is no hope 
to reconstruct a well-defined scalar-tensor theory. This 
may reflect the fact that the DGP model contains ghost 
modes around the self accelerating solution [g^] • 

Playing with the parameters 7 PPN and /3 PPN , we can 
reconstruct well-defined scalar-tensor theories for some 
values. Indeed, there is no insurance at this stage that it 
really describes a DGP model. Anyway, it would mean 
that a DGP model and the scalar-tensor could share the 
same background and perturbation evolution. Of course, 
this does not mean that DGP-5D are scalar-tensor theo- 
ries which would require a complete mapping of the de- 
grees of freedom of each theory. Also, the two theories 



cal value of which is determined by a Cavendish-like experiment 
today [see Eq. 1521 for the similar issue in scalar-tensor theory] . 
Without more information, we have to consider it as a pure bare 
parameter. The Poisson equation for DGP-5D takes the form 

A* = |n m0 ^-^(i + 4)- 

2 G cav0 V 3/3 J a 

with O m o = 8 7rG C avOPmo/3iJn . It follows that we can conjecture 
(see also Ref. that G C av = G(l + 1/3/3) so that 

This points out a difficulty. When we compare the background 
dynamics of DGP to another theory without having a full under- 
standing of the weak field limit, we cannot be sure that what we 
define as C m o in each theory refers to the same quantity. Fortu- 
nately this does not change the conclusions of § III. A. 2 but can 
affect those of this section. Note also that it was shown |8ljl that 
in the Solar system the gravitational constant was shifting from 
G to 4G/3 when moving away from the Sun. While the effect of 
DGP on light bending ans periheli precession have been worked 
out, the interpretation of the PPN parameters in this framewok 
is not settled yet. The same issue needs to be studied in any 
model of the classes C and D. 



can be hoped to be discriminated by other features that 
are not discussed here such as the high-rcdshift proper- 
tics (CMB, BBN), primordial predictions, strong field ef- 
fects (black holes and gravitational waves emission) and, 
as discussed above and illustrated on Fig. 111! local test 
of gravity. While all these predictions are available for 
scalar-tensor theories, this is not yet the case for DGP 
models (see table QJ. 

V. CONCLUSIONS 

In this article, we have discussed various physical mod- 
els of the dark energy sector that can be the cause of the 
acceleration of the universe, following a physically ori- 
ented classification relying on the nature and coupling of 
the degrees of freedom of dark energy as starting point. 
For each of the four classes, we have described specific 
signatures and tests (see Fig.^l. 

Beside the evolution of the background, we insisted 
on the central role of the perturbations to distinguish 
between the four classes of models. Restricting to the 
Newtonian limit, we propose to characterize them by a 
new set of functions besides the equation of state. Ac- 
cording to this classification scheme, the determination 
of the nature of dark energy can be outlined following 
a series of consistency checks in order to progressively 
exclude classes of models. Also, it can be used to quan- 
tify the departure from a pure ACDM in each class and 
to define target models in each class. Maybe it will just 
let us with the initial cosmological constant problem to 
face 0| after alternatives offered in the context of (well- 
defined) field theories have been exhausted. Let us recall 
that this analysis let a large class of solutions unexplored 
and also calls for tests of the Copernician principle. 

As exercices, we have constructed models from differ- 
ent classes that share the same background dynamics 
and eventually the same growth history of density 
fluctuations (see table [H] for a summary). These models 
may be differentiated on the basis of local experiments, 
strong fields effects and of large scale perturbations 
(including CMB). This illustrates the importance of the 
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theoretical prejudices when parameterizing the dark 
energy and the complementarity of data sets. 
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